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ABSTRACT 


Hallén type integral equations are formulated for an arbitrarily 
bent monopole above an infinite perfect ground plane and, by numerical 
solution techniques, current distributions are obtained. The coupled 
integral equations are founded upon the familiar electric scalar poten- 
tial and magnetic vector potential. Numerical results are presented 


for several cases including different bent angles and are compared where 


possible with available data. 
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{L. INTRODUCTION 


A linear antenna is one in which a "long dimension" can be discerned; 
for example, a long circular cylinder the length of which is much greater 
than the diameter is looked upon as a linear antenna. Initially, the em- 
phasis will be limited to the so-called cylindrical antenna driven at its 
mid-point by a voltage source V applied to a vanishingly short gap at the 
antenna center. The antenna will be considered to reside in a lossless , 
homogeneous and isotropic medium of infinite extent. | 

A thin wire antenna is commonly defined as a linear antenna which has 
a radius much less than a wavelength. There are several approximations 
which may be employed in the analysis of thin wire antennas. 

(1) The current on the cylindrical surface is assumed to be 
entirely in the axial direction. 

(2) The current distribution is assumed to be constant around 
the periphery of the wire at a given cross-section and hence can be 
approximated as a filament of current on the axis of the cylinder. 

(3) The boundary conditions are applied only to those components 
of the total electric field which are tangent to the wire surface and 
Paeaiiel to the wire axis. 

Many investigators have considered the theoretical problem of deter- 
mining the current induced on a thin wire scattering antenna. [Ref. 1,2,3]. 
These investigations have been extended to include the treatment of ar- 
bitrary configurations of thin wire scattering antennas intersecting at 
a single junction. [Ref. 4]. Three particular configurations for which 


numerical results have been obtained are a set of nonintersecting thin 





straight wires, [Ref. 5], a set of perpendicular intersecting thin straight 
wires, [Ref. 6], and a pair of skewed crossed wires. [Ref. 7]. 

One general approach to finding the current on thin wire antennas is 
to obtain a set of integral equations which involve the unknown current 
distributions. Using the equation of continuity, one may determine ex- 
pressions for the vector and scalar potentials in terms of the current 
distribution. It is in obtaining these integral equations that a choice 
can be made. 

Following the Pockington equation formulation, the expressions for 
the vector and scalar potentials are used along with Maxwell's equations 
to derive a system of integro-differential equations for the unknown 
current distribution. An auxiliary scalar function is defined and intro- 
duced into the system of integro-differential equations. The solution of 
these now ordinary differential equations subject to appropriate boundary 
conditions gives the desired set of’ integral equations for the unknown 
@emrent distribution. 

If the Hallén integral equation approach is used, one independently 
establishes two expressions for a given component of the vector potential 
in the direction of a wire element and these expressions are evaluated on 
the wire surface. One expression is obtained from the solution of an 
inhomogeneous differential equation satisfied by the component of vector 
petential. 

Another expression for vector potential is obtained directly from the 
familiar vector potential itself. These two expressions must be equal 
and, hence, are equated. Next, the boundary condition that the tangential 
component of electric field on the surface of the wire element must be 


zero is applied. Now the arbitrary constants of integration which arise 





in the solution of the inhomogeneous differential equation are evaluated. 
The constants of integration must have assumed values which are in com- 
memance With all other required conditions. [Ref. 8]. 

There have been numerous methods proposed for solving the integral 
equations, however obtained. These methods have fallen into a few broad 
catagories: 

Gp Srteratrons (Ret. 9. 10]. 

(2) Expansion of the current in a Fourier series [Ref. 11]. 

(3) fhe variational method [Ref. 12,13,14]. 

(4) The Wiener-Hopf techniques for the case of the long 
emeenna [Ref. 15]. 

(5) Numerical methods [Ref 1, 16]. 

One effective numerical method of solving the integral equation is 
the so-called direct integration technique (Method of Moments). [Ref. 16, 
17). Basically the technique utilizes an expansion of the unknown func- 
tion in some convenient manner. This expression may be a finite sum of 
piecewise continuous functions which are non-zero over some small inter- 
val of definition and vanish elsewhere. A substitution of this expres- 
Sion for the unknown function into the system of integral equations 
reduces them to a system of linear algebraic equations for the unknown 
coefficients. These linear algebraic equations may be solved by conven- 
tional techniques. 

The purpose of this thesis is to formulate the coupled integral 
equations for an arbitrarily bent monopole, above an infinite, perfectly 
conducting ground plane, in a Hallén-type integral which is founded 
entirely upon the familiar scalar potential and magnetic vector potential. 


These integral equations are then solved using the direct integration 





fechnigue. Throughout the thesis it is found that the coupled integral 
equations obtained by the Hallén approach can be readily obtained with- 
out employing an auxiliary scalar potential. How the various arbitrary 
constants, which arise in the development, are related to the electric 


scalar potential is clearly shown. 
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II. DERIVATION OF MATRIX EQUATION 


In this section the system of integral equations for the current 
distribution on an arbitrarily bent monopole, above an infinite, perfect- 
ly conducting ground plane, driven at its base by a delta gap generator 
is derived. From the equations, an operator is defined such that the 
method of moments can be applied to find an approximate representation 
of the operator by a matrix. The solution for the desired current dis- 


tribution is obtained by inverting the matrix. 


A. FORMULATION OF THE PROBLEM 

Figure 1 defines the coordinate system. The antenna consists of two 
thin wire segments with a single junction. The segments have radii a, 
and a. and lengths hy and he: The primary source is the field in the 
gap which is maintained at V/A volts per meter by the generator. The 
induced source is the current on the dipole. The solution is carried 
out aSSuming the wires are perfect conductors. With this assumption the 
tangential component of the total electric field at the surface of each 
wire segment muSt vanish. Due to axial symmetry of each wire segment, the 
phi-component of each electric field must vanish. Hence for each wire 


segment the following is obtained: 


(1) 


ll 
ee 


ly 8) ae 
EO (ry) + EY (r.) 


(2) 


ll 
© 


Ls pe 
Ba'e,? * E,@,) 


where EO (r.) is the tangential component of the induced electric field 


at the surface of the z-directed wire and EF (x) is the tangential 
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Figure 1. Antenna Geometry . 
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component of the primary electric field at the surface of the z-directed 


wire. Equation (2) is similarly obtained for the s-directed wire. 


B. DERIVATION OF THE COUPLED INTEGRAL EQUATIONS 

The Hallén equations are derived from a classical viewpoint for the 
case of a cylindrical antenna driven at its center by a so-called slice 
generator. [Ref. 9, 18, 19]. The general procedure for developing a 
Hallen type equation is as follows: 

(1) Determine the magnetic vector potential component from the 
potential integral equation for A. 

(2) Determine the magnetic vector potential component along the 
cylinder from the inhomogeneous differential equation which the magnetic 
vector potential component satisfies. 

(3) Write the desired integral equation by equating the two ex- 
pressions obtained in steps one and two. 

(4) Apply appropriate boundary conditions on the surface of the 
cylinder to the equation obtained in step three. 

1. Vector Potentials 
Consider the antenna shown in figure 1. Assuming harmonic time 


jut é ; 
dependence (ed ), the vector potential due to the antenna current is 


as call -jk|r-r'| 2 
Ge) 322) 3) eT Cap (3) 
; At |r-x' 
ae r-r 
«ee be ee 
CP maaan Je dV (0' ) (4) 
At a \p~ 7 


Due to the geometry the vector potential has only a z-directed and s- 


directed component which can be written on the axis of each wire segment 


Ale 





as 








r a7 JKR 
a. pee ' 
A, (2) in x, ee) dz (5) 
2! 
r ike 
= —— y ! 
A_(s) an | I (s ) ds (6) 
ea! 
where R= |r-r'| and P = jJp-p'{. R and P were evaluated at r' and 9' 


equal to the radius of the respective wire segments with ¢' equal to zero. 
2. Non-Homogeneous Differential Equations 
It is now seen that A, (z) and A (s) can not be determined di- 
rectly from equations (5) and (6) without knowledge of I(z') and I(s'). 
A differential equation is now derived which can be solved for A, (2) and 
A. (s) 


— 


E = -jwA + 





jw erad (div A). (7) 


mewelectric fields produced by the currents are, according to equation 


(7), 


tt 





meee) jwA (E)U + - : erad faiv(a, Gu, )] (8) 


Zz Zz Z jwue 
(4,0, )]. ©) 


Having applied the gradient operation in circular cylindrical coordi- 


I! 


E_ (p) 





nates and with axial symmetry eliminating the phi-~component of E “(r) 
and E. (), the fact that the total tangential electric field on the 
surface of the antenna must vanish is applied. The following results 


are obtained: 


2 
er G. i. A, | (10) 


Il 


2 
5 = age 
<7 [A &) | #k A) 


2 
jk BOO | On 0 
Ww Bg (P .) Os E a @,) |» oe) 


2 — —_ 
= [4,@,)| +k A. (® .) 
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where =. and p. denote any point on the z and s-directed wires respec- 
tively. Equations (10) and (11) are two coupled partial differential 
equations because of the last term on the right in each equation. A 
Hallen solution of these two coupled equations is needed. 


pee OluLttlon sor Coupled Dirterential Equations 


Equations (10) and (11) are of the form 


2 
oie) $ ten = u(x) + oO) ; (12) 
dx 


The general solution to equation (12) was contained in Ref. [235 


was found to be 


X : 
A < a rd 
mic) = A cos kx + B sin kx + an Ke ) + uae sin k(x-x )dx . 
“4 (13) 


The integral portion of equation (13) is zero at x equal q and its deriv~ 


aelve is zero at this position. One integration by parts results in 


x 7 
i RG? sin k(x-x ax’ = -V(q) sin k(x-q) 
eC 


= 
+ k il V(x’) cos k(x-x’)dx’ Gks) 
/ 
=a] 
which, when put into equation (13), gives 


f(x) = A cos kx + B sin kx - — v(q) sin k(x-q) 


ie = a f 
a L U(x’) sin k(x-x‘)dx’ + V(x’) cos k(x-x’)dx’. 
a i 
—q =q 
(15) 
Since V(q) is fixed, the first three terms are combined and new constants 


Oeeained are 


a= At< V(q) sin kq, and (16a) 


b 


I 
bd 
| 


= v(q) cos kq. (16b) 


iS 





Combining equation (16) with equation (15) gives 


fiejye= a cos kx + b sin kx + * | U(x’) sin k(x-x’)dx’ 
ic 
s 
+. L V(x’) cos (x-x’)dx’. Cle) 
—q 


With equation (17) available the solution to equations (10) and (11) is 


‘determined with 


ZL 
fe) — jk —- ; 
as A, (tr) = V(x); oe P(r) = U(x) for equation (10) 
a2 
and ao A (0_) = V(x); ag Pip) = U(x) for equation (11) 
dw ws ? W Ss S q 


which gives 


: Z 
A (z) = a_ cos kz # b~ sin kz Bs I E (mn) sin k(z-n)dn 
Z Z Z | WW 7 Z, 
= 7, 
fe) 
Z 
oer ie = A (n) | cos k(z-n)dn and (18a) 
6) 


S 
A (s) = a_ cos ks + b_ sin ks + x | E (mn) sin k(s-n)dn 
S S S WW) ete 
0 


S 
- —- E A (n) | cos k(s-n)dn. (18b) 


The Sa] means A, (x) is differentiated with respect to y and then 
has to be evaluated along a z-directed line on which it is dependent. 
This is similarly true for = A, (s) | The arbitrary constants of equa- 
tron (18a) are ai» b. and Zo and any one of these three may be specified 
directly, but the other two must comply with physical boundary condi- 
ditions. This applies to equation (18b) and constants ao: b. and S 


also. 
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tne ocalar Potential 
With equations (18a) and (18b) one faces the problem of what to 
wemoveut the constants of integration. One can enforce the boundary 
condition that the scalar potential must be continuous at each wire seg- 
ment junction to eliminate a constant or relate two constants. To obtain 
such an expresSion is the next task. 
Differentiating» equation (18a) with respect to z and (18b) with 


respect to s according to Leibnitz's rule [Ref. 24] gives 


Z 
A (a ookanesin kz + kb cos kz + cls Hein) sin k(z-n) dn} 
OZ Zz Z, OZ W) Zz 
n=Z 
O 
eee Ay (z) +k [2 A (n)] sin k(z-n)dn (19a) 
oy vey 
and 
S 
a A (s) = -ka_ sin ks + kb_ cos ks + mo els E (n) sin k(s-n) dn} 
os s S S oS W) S 
n=s 
O 
3 
3 ‘tO 
> A _(s) qe [2 a, oo) sin k(s—-n)dn. (19b) 
n=s 
O 


fhiesLorentz condition relates the magnetic vector potential A to the 


electric scalar potential and yields 


= ic, @ = 2a @+—2a@ d AO 
me yy dz 2 ss a 
Z 
_ J ma 
pre) = a. A _(o) 1 ere date (p) . (20b) 


Rearrangement of equations (19a) and (19b) gives 


Z 
ee : 
- pet 2 (5) = -ka sin kz + kb cos kz fee Sus [r (n) sin k(z-n) dn} 
W z Zz OZ W Z 
n=z 
Zz O 
,) 
[— A (n)] sin k(z-n)dn (21a) 
x {2 a] 
n=Z 


ay 





and 


S 
2 
mak .s a: Soak fe : 7 
a (s) ka. sin ks + kb cos ks + 56 { 7 [za sin k(s n)dn} 
n=s 
O 
S 
e : 
+k [so 4, (0) sin k(z-—n)dn. (21b) 
n=s 
fe) 


which is valid along the z- and s-directed lines respectively. 

6” (z) Pee vesuelalcealar potential alongethe z—- directed wire 
and $° (s) is the total scalar potential along the s-directed wire, both 
of which are due to all charges on the wire segments, but the potentials 
are expressed in terms of the constants of integration. The partial 
derivatives expressed in equations (2la) and (21b) must now be evaluated. 

Going back to equation (13) it is seen that the integral portion 
of (13) is zero at x equal q and its derivative is zero at x equal q. 


Now in the derivation these forms are used 
x 
> aa 


Their relationship to the antenna problem is given below equation (17). 


x 


U(x?) sin k(x-x*)dx” + { Vix) cos Kix-x )dx . (22) 


ato 


( 
=q xsd 


Now applying these facts the result is 


Z 
“73 jk ; = : | 
A { - | EY {n) Sik CZ n) dn} = 0 for all Zo (23a) 
n=Z 
O 
Ss 
3 7k z 
= { ‘ | E_(n) Sani k(s-n)dn} =) kee Gull! s° (23b) 
n=s 


This result implies that one may choose s_ and Zz, as desired, but at the 
O 

Same time doing this in a manner to lessen subsequent calculations. Once 

Sy and z, are determined the other two constants of integration in each 


equation must meet boundary conditions. 
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C. BOUNDARY CONDITIONS 
The application of these necessary boundary conditions to the cur- 
rents on the thin wires should provide the additional equations to yield 
a unique solution for the current distributions. 
ie At The Junction 
As has been stated by Taylor [Ref. 4], there are essentially 
two conditions which must prevail at all junctions: 
(1) Kirehoff*s current law must hold. 
(2) Continuity of the scalar potential must exist. 
The first of these requires that there be no buildup or accumulation of 
charge at the junction. Stated differently, the total charge entering 
the junction must be equal to the total charge leaving the junction. 
The second condition requires that the potential on each of the wires at 
the junction be the same. There can be no jump or discontinuity in sca- 
femepotential at a wire junction. 
2. At The Free Ends 
According to previous considerations of thin wire theory [Ref. l, 
2,4], the current must vanish at the free ends of all wires in the con- 


figuration. 


Ihe 





ieee er reCAT IONS ARISING FROM GEOMETRICAL SYMMETRY 


Any symmetries in the geometry of the antenna elements are reflected 
in the structure of the generalized impedance matrix fee . Judicious 
use of these symmetries can result in considerable reduction in both 


computing time and the number of memory locations required for a solution. 


A. IMAGE SYMMETRY 

Referring to figure 1 the antenna has image symmetry about the plane 
z = O through the midsection of the antenna. The boundary condition at 
a perfect electric conductor is zero tangential components of electric 
field. An element of source plus an "image" element of source, radiating 
in free space, produces zero tangential electric field over the plane 
through the midsection of the antenna. According to uniqueness concepts, 
the solution to this problem is fle the solution for a current element 
adjacent to a plane conductor. The necessary orientation and excitation 
of image elements is summarized in many texts. The main point though 
is that symmetry must be maintained in the image problem. Once image 
theory is used and the resulting symmetry is recognized, all that is 
needed is to find the kernel for each "image" current element and com- 
bine it algebraically with its corresponding "“non-image" current element 


kernel which reduces the number of required equations by a factor of two. 


freee (ORPLITZ PROPERTY 
The Toeplitz property of the coe matrix can be used to realize 
considerable economy of computing time. An element of the matrix depends 


only on the absolute value of the distance between the source point m 


20 





and the field point n if there is no bend involved. Since there are 
only N distinct distances involved, only N different quantities need 
be calculated for a matrix containing N x N elements. These N quanti- 
ties are calculated first and the computer is instructed to arrange 
them in the proper locations to form the complete matrix. 

The antenna has a bend so only certain portions of the matrix can 
be obtained using this property. It is the image of the bend that pre- 


vents the entire matrix from exhibiting this property. 


iad) 





IV. NUMERICAL SOLUTION 


One of the most general and most useful techniques for solving pro- 
blems in antenna theory is the direct integration technique. A detailed 
discussion of the technique has been presented by Harrington [Ref. 16]. 
Since this basic technique is employed throughout the remainder of this 
thesis, a short quantitative discussion using the notation of Harrington 


is included at this time. 


A. METHOD OF MOMENTS 
The equation to be solved is of the form 
CE) =) 2 (24) 
where L is a linear operator, g is the tangential component of the applied 
field on each conductor and f is the unknown current to be determined. 
Using the terminology of desaeon nee 17] it is evident from the 
uniqueness theorem that the problem is deterministic. Since L and g are 
given, it is a problem of analysis. 
A unifying concept for the solution of such problems is afforded 
by the method of moments. In order to apply the method, it is necessary 
to define: 
Gl A set of expansion or basis functions fi> Eayeees fT in the 
domain of L. These known functions must be linearly independent and 
chosen such that a reasonable approximation to f can be obtained by a 


finite series expansion of the form 


N 
‘oy ) aa (25) 
n=1 


a 





where the as are unknown constants to be determined. If equation (25) 
is substituted into equation (24), one obtains 


N 


bth Ey) =p. (26) 


(2) A set of weighting or testing functions Was Woos--oW , Wik? 


...,W,, in the range of L and linearly independent must be chosen. If 


N 
the same set of functions is used for expansion and testing, the pro- 
cedure is known as Galerkin's method and is equivalent to the Rayleigh- 
Ritz variational procedure [Ref. 17, 25]. 


(3) An inner product suitable for the problem. For most electro- 


Magnetic field problems, a suitable inner product is 
2 
<W,f> = | f(n)W(n) dn! (27) 
a 


where ny and n, are ElescimOven TENG region Of Gmeerest. The resulting 


expression is 


i et ba 


in Ce = We hoe ee eee Ns C25) 


n 
Inspection of equation (28) shows that it can be written as a matrix 


equation by defining the following matrices 


eee 4. Wes 


Ie ik tea. 1 nN 
< Is «#0 @ 
ere Wo » i <W,,Lf,> <W,> Lf. ons 
mn : Y . 
“Wyo bE > Wyo LE, > pate <Wyy> LE? 
? 7 la] = Oy , and an = <W, 58> : (30) 
oo) <W, 8? 
Chay <Wyy 8? 


Z3 





The matrix equation may now be written as 


Bell me : ue : (31) 


By analogy with multiport network theory, (ees is called a generalized 
impedance matrix and Md a generalized voltage matrix. Only if [za 
is nonsingular does its inverse exist and the desired solution is 


a = 
fo} = (2 1° (v1. (32) 


Since the a's are then known, the solution for the unknown current f 
is determined by equation (25). 

The expediency of obtaining a solution for f depends not only on L 
and g, which are uSually determined by the physical problem, but also 
on the functions a and Wee Extreme care Should be taken in choosing 
these two sets of functions. Experience with problems of the general 
type considered in this thesis suggest that only certain particular 
choices for the expansion functions and for the testing functions yield 
equations which are amenable to solution. In particular, when the test- 
ing functions are chosen to be Dirac delta functions, the resulting 


technique is usually referred to as point matching. 


B. EXPANSION FUNCTIONS FOR THE UNKNOWN CURRENT 

There are several popular choices for the expansion functions which 
have been considered in the literature. Some of the more common of these 
are piecewise constant, piecewise linear, piecewise parabolic and piece- 
wise sinusoidal [Refs. 20, 21]. All of the aforementioned fall under the 
general label of the method of subsections. This method has the antenna 
contour divided into a number of intervals or subsections where the 


boundary condition of equation (1) and (2) is to be applied. The lengths 


Ze 





of the intervals need not all be equal. For certain types of problems 

it could be advantageous to weight the lengths of the intervals. In order 
to keep the formulation as simple as possible, however, the present 
solution uses N-]l intervals of equal lengths and two intervals of half 

the length of the N-1 intervals. These half subsections were to allow 

a match point to be placed at the ends of the antenna and one on the 

bend of the antenna. In addition to these piecewise functions is another 
family of expansion functions which are commonly called entire-domain 
expansions. A discussion of the use of these entire-domain functions 

will not be attempted in this thesis. 

All of the calculations and results presented in this thesis are 
obtained by the method of subsections in conjunction with point match- 
ing. In particular a series of piecewise constant and piecewise sinu- 
Soids are used to represent the unknown currents. The resulting 
equations are forced to be satisfied at a set of discrete points called 
match points. Using Dirac delta functions as weighting functions ef- 
fectively requires that equation (24) be satisfied at N points on the 
contour of the antenna. 

Harrington and Mautz [Ref. 2] have solved the single straight wire 
problem uSing three different procedures, all within the scope of the 
method of moments. In order of increasing computational complexity, 
these include: 

(1) Point matching with pulses as basis functions 

(2) Point matching with piecewise linear basis functions 

(3) Galerkin's procedure with piecewise linear basis functions. 
As the number of subsections increases, the second solution converges 


about twice as fast as the first. With antenna segments less than A/10 
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in length, no significant difference in results is observed between 
Galerkin's procedure and the second procedure. Using piecewise sinu- 
soids is equivalent to the second method in speed of convergance and 
and accuracy, but piecewise sinusoids do not require the segment length 
to be less than 4/10 in length to obtain this accuracy and speed. 


Using Dirac delta functions as weighting functions gives 


6(z - z_)U, (33a) 


TM 


W 
S 


I 


é6(s - s_ JU. (33b) 


dihew ra ch Je is a unit vector directed along the positive z-axis and U, 
is a unit vector directed along the positive s-axis. Using the piece- 


wise constant subsection gives the following expressions for current 
N A 
Gam) IP (z)U (34a) 
nn Z 
n=1 
I(s) = ; ies) ie (34b) 
q=1 qq Ss 
where I and He are the complex constants to be determined and the pulse 
functions are constant over the increments. This statement in equation 


form is 


1 ZuaZ az 
n 


pace ntl 36%) 
= 0 elsewhere 

1 Ss <s<s 
p(s) = q qtl (35b) 


0 elsewhere 
Similarly, the piecewise sinusoid expression gave the following current 
equations 


‘ t of 
N I sin k(z Zee Len 


Sinek(g ez) 
Ce y ntl n 
n=l 


J Caine (Slee 
h Zz 


‘ Fis ol EY 
sin k(z hay z_) 
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eectoeloas'  —s') +t sin k(s'~s') 
Sq qt 


+h 1 


Q 
Hp) 


; i Moe (36b) 


1 e 1 
sin eens 


“sD 
C. LOCATION OF MATCH POINTS 

There has been little quantitative study done on the precise effect 
of the location of the match points on the solution; however, it has 
been the experience of this author that the location of these match 
points significantly affects the accuracy and convergence of the solution. 
The number of match points and their location depends on such factors as 
the form of the integral equations (Pockington or Hallén), the type of 
current expansion used (pulse, sinusoids) and the allowable time and 
core for computation. 

Extensive calculations have been carried out [Ref. 22] where the 
match points are located at onde of each current subsection. The results 
of these calculations indicate the solutions are quite stable and reason- 
ably curate even for problems where a large number of match points is 
used. 

Other investigations have indicated that a good choice for the loca- 
tion of the match points is at the center of a current subsection. As 
the author Migeovered these conclusions have to be tempered with other 
information. One question is how and what kind of experimental verifi- 
cation is used. If a scattering cross section is used, the demonstrated 
results may differ from those obtained for an excited antenna. The reason 
is the cross section calculations effectively average the source and 
current Peet ieas Si the structure. Thus, a convergence study 


pertaining to input admittance, using scattered results, may be somewhat 


ey 





inconclusive since the segmentation and match point dependence are 
directly related to source model and numerical convergence. 

All of the results presented in this thesis are the result of locating 
the match points at the center of a current subsection for piecewise con- 
stant expansion and at the ends of the current subsections for piecewise 


sinusoid expansion. 
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V. APPLICATION OF THIN WIRE THEORY TO AN 


ARBITRARILY BENT WIRE MONOPOLE 


A. GENERAL CONSIDERATIONS 

feeeneral theory for determining the currents induced on an arbitrar- 
ily bent wire monopole has been outlined. The purpose of this section is 
to illustrate how to obtain the exact solution. The configuration to be 
studied is composed of two straight nonorthogonal elements intersecting 
at a single junction which are driven above a perfect infinite ground 
plane. The problem is developed in a manner which allows general var- 
iations in many of the fundamental parameters. For example the lengths 
of the individual elements may vary, the angle between the individual 
elements may vary and the radii of the elements may vary. It should be 
noted that all variations in parameters are still subject to the basic 
assumptions and restrictions imposed by thin wire theory. These restric- 
tions do not greatly inhibit the usefulness of the investigation since 


much general knowledge can be obtained from the thin wire approach. 


B. COUPLED INTEGRAL EQUATIONS 
The coupled integral equations for magnetic potential which were 


developed in a previous section are 


Z 
heez w= 9a. cos kz + b. sin kz Fs | E (n) sin k(z-n)dn 
Z Z vA W Zz 


Zz sal 
-i| [a A) cos k(z-n)dn (37a) 


n=Z 
oO 
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S 


A (s) = a cos ks + b~ sin ks + ik | E (n) sin k(s—n)dn 
Ss Ss Ss Ww Ss 


Ss as 
1 3 
-< L=-[A_ (m2) }} cos k(s-n)dn (37b) 
n=s 
O 
and the scalar potential equations are 
D Z 
mee (5) = . | 22 a ae 
a » (Zz) ka, Sin kz + kb | COs "kai | ty lA @ Sin k(z-n)dn (38a) 
n=Z 
and 
Pe Z 3 
an @ (Zz) = ~ka, sin ks + kb. cos ks ek i — Sin k(s—-n)dn. (38b) 
n= 


The constants a. bo a» and D. must assume values which render equa- 
tions (37) and (38) in compliance with the following boundary conditions: 

(1) The scalar potential,$, must be continuous along the wire 
segments except at the delta gap generator and in particular the scalar 
potential must be continuous at the wire junction. 

(2) The magnetic potentials, A. and Ai. must be continuous along 
the wire segments and at the wire junction. 

(3) Current must vanish at the free end of the wire. 

(4) The axial currents must satisfy Kirchoff's current law at 


the wire junction. 


The magnetic potentials, A. and An developed earlier, are, 





kR 
A(z) = ah I (z' yea (39a) 
Aa 
i Bat 
and NG) we | Ij(s") ds‘ (39b) 
er 
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where R= |r - r'| and P = |p - 0'|. R and P were evaluated at r' and 
p' equal to the radius of the respective wire segments and $' equal to 
zero. Using these values for |r - Ta and lp — o'| and assuming the sur- 
Face current flows as an equivalent line current along the s and z axes 
is equivalent to using the reduced kernel. This will be proven in Ap- 
pendix A along with conditions when the approximation is no longer valid. 

Looking back at the development of equation (37) it is seen that the 
mathematical formulation was given without any physical insight as to 
where in space the equation was valid. In going from equations (10) 
and (11) to equation (18) there is a built-in assumption not expressly 
stated. This assumption is that the radius vector is fixed at some value 
which results in equation (37) being valid only along a line parallel to 
z ands directed wires. This result will be fully shown in Appendix B. 
Since equations (37) are valid along any line parallel to the z ands 
directed segments, one can require ee to be along the surface of the 
wire segments which is along a line (a,o,z) and (a,o,s) respectively. 
This surface is also where the tangential E field is zero. Since the 
Source points and field points could not lie along the same line, the 
reduced kernel has to be used with its equivalent line current flowing 
along the axis of the two wire segments. 

Figure 1] de Bae the co-ordinate system which is composed of two 
co-ordinate subsystems. They are the (X,Y,Z) and (P,W,S) co-ordinate 
systems. The kernels for the wire segments are then expressed in these 


co-ordinate systems. 


a7 3kR” ; -4kR 
K (x,y,z,z') = ——— ka ec, 2,2.) = (40a) 
Z Z Z i Z 
R Re 
ae 
~jkR° F -jkR; 
K (x,y,z,8') == Kecx.yoz>s'!) = (40b) 
S S S i S 
R Re 


Sil 





eS Ss Soi 


A 


Ue 


{1 


{{ 


if 


ft 


A A 
= ~U, Cos20+ ue SiN 20 


A A 
-UL Sine + Uy Cos © 


A 
Uh co. 8 +Uy Sin @ 
A A 
U, SiIN@ +t U, Cos © 


A A 
U, Cos 8 — Ue Sin 0 





Figure 2. Antenna Coordinate Systems 
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eZ cee 
a Jkr ; eo jkr; 
! = ! Ss fe See 
G (p,w,S,2') = = G(p,W,S,z.) - (40c) 
ie ie 
it 
ov dk of ov ik 
' = ' Sa 
G_(p,W,S,5 ) = - G(p,W,5,5,) = : (40d) 
18 Lr. 
ah 
i 
Pies ey = ye Ge = 2)" 
i 
i Z 
r= {(p -p') + Ww- mye (Ss = Ss ) i , in general. 


In order to account for the images due to a perfect ground plane one 
must include the image current kernels which are denoted above by the 
Superscript 1. 


With reference to figure 2, 


hy 0 
7A. ee + 7 ! ! oye 12 ! ul ‘yds! 
A(r) v4] _ (2 )K (x,y, 2,2 ddz' 4 | 7 b2,IK, 22.25) zt 
= Te 
he gar oi a 
‘ '\yv t f - : ! aa eee ! ' 
+ Un | Ii(s )K(K,y,258 )ds Ue; | To (s,)Ko(s,y,2,8,) ds, (41) 
s'=(0 si =-h 
a S 
and hb, 0 
Nees i ' ' ' eee - ' ' 
A(o) v, | I (z )G (PsW,S,z diz it [ 1,@)6,@.¥,s,2))42/] 
z'=0 Zeal 
a ne A : i 
' ' pas ' ' ' y 
+ U, I_(s')G. (p,w,s,s')ds Ue. | I (s,)G.(p,w,s,s;)ds.. (42) 
s'=Q =n 


al S 

When the Uns and UL unit vectors are broken into their respective U, 
and U, components, they give the required A(x) and i Similarly 
breaking U, and Uns into their respective U. and U. components gives 


A, (°) and A (0). These four magnetic potentials are needed in equation 


37. These four required magnetic potentials are: 
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Z = 
BA) - | TC Gay)z, 2") dz’ + cos6 | I. (s')k° CR ihe ZS as) » (43a) 
z'=0 cone 
h 
An = ° 1\L) ' ' 
a se) = sin @ I (s WitiGey. 25° )ds", (43b) 
s'=0 
as 
An 1 
7 ee zs | I_(s'){G_ (p,w,s,s")ds! COS cd G_(p,w,s,s')} ds’ + 
s'=0 
h 
Z 
cos 8 | I (2! Nex Cie wise, )dZ. (43c) 
¥— 
and ‘ : 
At . 1 ' ' 
ee <= —_— 1 ' — 
i A (Pp) Sade 20 | I (s')6 (p,w,s,8')ds 
s'=0 
h 
zZ, 
Sim 0 | (2! NG (paWweerme Jaz, (43d) 
z'=0 
Placing equations (43) into equations (37) gives 
h h. 
z S 
i jt I (2')K? Gaye 2,2 )dz + cos 0 {x a (Se Ke (x3y52,5 ds’ = a cos kZ 
“a s'=0 
Z, z, 
b sin kz Ai gpils E (n)sin k(z-n)dn - 4 eae (n)}} cos k(z-n)dn (44a) 
Z Ww Z K dy y 
n=zZ N=Z 
O O 


Z 
I @ "CG (.w,8,2' dz 


——— 
= 


a 
S 


S Z 
Vs =” z'=0 
a_cos ks + b sin ks + ake E (n)sin k(s—-n)dn - | {*[A,(n) }} cos k(s-n) dn. 
n= 


(44b) 


Examination of equations (44) indicated that the partial derivatives 
of A, (r) and A Co) are needed. The results for HGS) Wilieepe shown and 


A Co) Wittivrollow ine a) Similar manner. 
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ihe form of a is 





on SKR 
A. (r) = 5 (454) 
dA (r) eC Ochs 
Bs y Sy et) Salleh og 
oy oR oy 
3A —jk 
Ae NEOCIESER Gy ye (45b) 
ay R? 
Similarly, for A); 
PAO) + jkr) (w= whe IK sie? 
ow ae 


Placing equations (45) and (46) into equations (44) and examining 
the results closely, it is seen that the double integrals have an inner 


integral of the form 


cee “ees cos k(8-n)dn (47) 
2 3 
R R 
i 
where R = (ee + g(a) + [n - f(a) )*}° 
and f and g are any function of the parameter a. Reference [23] gave 


the exact solution for this integral which reduced the double integral 
to a single integral and thereby greatly reducing the computer inte- 

gration time and numerical errors. The exact derivation will be shown 
in Appendix C. Now place the result of Ref. [23] into equation (44a). 


With rearrangement of equation (44a) the following equation is obtained: 
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Eile, cos kz + [23] we sin kz + (A) | E (n) sin k(z - n)dn - 
U Z Ul Z Uw Z 


n=zZ 
O 
Z, 
“ [F. en (x? Ga ake 15) (LEO eee, - ie tana Jas' 
2 Ce iS) )) | 9) — 
O 
Z 

| sx - cos 6 (i (s' )K° Ceayac.s jds™ = 0 (48) 
Z'=Z s'=s 


femilarly, for equation (44b): 


S 


An An : An 4k , 7 
[ a a. Cos ks + [ b. sin ks + irra | E.(n) sin k(s - n)dn + 








n=s 
he i S 
sain 26 ex 5") (Gaim eS ) ]cos = n) - = k(s- al) )ds’ an 
ee r(x* So a(S) a =) wee 
h ° See 
; Ss 
ae | 1.(s') (6”(p,w,n,2') {tote Hees Reson) — jrsin ESM)! Jast - 
ee (fe SrA) 2 ae 
h 2 i 
Z il 
| I. (s') [G_(p,w,s,s') - cos G_(p,w,s,s') ]ds' ~ 
s'=s 
O 
Z 
cos 6 | I (z")G°(p,w,s,2")dz' =Q . (49) 
ZtZ 


For a Structure driven by a delta gap source of emf, V, on an element 
at a location displaced g units above a ground plane the following equa- 


tions are obtained [Ref. 20}: 


V_ (2) - [- ~Vé(n-g) sin k(z-n)dn (50a) 
n= 
Vf) = Bey sin k(z-g) U(z-g)]. (50b) 
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It is desirable to obtain a representation of the emf when the 
Structure is base driven, i.e., when g is zero. A Structure above a 
ground plane, regardless of the nature of its excitation, always gave 
rise to zero scalar potential everywhere on the ground plane. This is 
a direct result of the uniqueness theorem and boundary conditions at a 
perfect conducting surface. Thus, it is necessary that, if the element 
and its source are to be equivalent, the impressed emf must be sym- 
metric; otherwise, the currents and charges in a free-space "equivalent" 
of the structure above a ground would not satisfy imagy theory require- 
ments. 

It is from this discussion that the modification to equation (50b) 
is found. When a base-driven structure is being analyzed, one sets g 
equal to zero and V equal to V/2. If V/2 were not used in place of V, 
one would be treating an antenna driven by 2V when the free-space 


"equivalent" structure is formed. These results are 


V (z) = - elie Sin kz. (51) 
Z Ww 2 


Since the s-directed segment is not driven by a separate generator, this 


gives 
S 


Vv. (s) = 0 or ik VOlecsanek(s = njdn = 0 . (52) 


_ 


U) 
= S 


O 

Now the boundary conditions are applied. Due to symmetry of the 

charges on the structure and its image, the scalar potential must be 

zero at the point where the vertical element enters the ground plane. 
This implies 


=o (0) = 0 (53a) 


which gives be sao (53b) 


oF 





fi 


Next by condition one the scalar potential must be continuous at the 


wire junction which gives 


°(h) = ¢° (0) (54a) 
and results in 
Zz 
-ka sin kh - kb +k Cas (n) ]}sin k(z-n)dn - et cos kh = 0 (54b) 
Zz Zz S oy y w 2 Zz 
n=Z 
O 


Conditions three and four require 


Oe and (55) 
Q 
Zz  -§ 


Clearing all constants and applying all boundary conditions, the fol- 


lowing equations for I) and I(s) are obtained: 








h Z 
S : ae 
oe Ree eee | I Ge Gay me) ( SS )Jecos k(z-n) - jRsin eigen) Jas’ 
ae : R(R” = [n-£(s')]*) 
s'=0 =0) 
h h 
iM t\y7 ' ' 2 t\yo ' ' 
- aa Sin kz — I DRG, y Zee )dz cos 4 I(s Wieetx.y,2,5 ds 
z‘'=0 s'=0 
= 0 | 7 (57a) 
he S 
[“Ta_ cos tage | 1G) © Gangs | EC }as' 
. Snes r(r° = [n-£(s')]*) - 
bh S 
ee 1p, Bee ey | I (2') (67 (p,wyn,2') {Ee Ea Lig k(s-n) - ins k(s- wy Ae 
> a - [n-£(2')]*) 
z'=0 [n- n=0 
S me 
- | I_(s')G_(p,w,s,s')ds! - cos 6 ese 
eo f= s'=0 
Zz 
- cos 6 i I (z')G’(p,w,s,z')dz' = 0 (5 7b) 


g'= 
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and 


h h 
= ' ee Z 
: AT a sin kh + | 1 fe" CS Gag. 2:90) eee aD ast 
5'=0 R(R’ = [n-£(s") ]*) =0 
i> - 2 cos kh. = 0 (57c) 
a6 s 60 Z 


= 0 (57d) 


ne oo (57e) 


There are five unknowns aa, Dane aa: I (2*), and I(s') to be found 
S 
and five equations which relate these quantities; hence, there is a unique 
solution to this system of equations for the five unknowns. [Ref. 20]. 


Obtaining the current distributions, all the other desired quantities 


een as €lectric field, input impedance can be determined. 


fee NUMERICAL SOLULION 

The system of integral equations given in the previous section must 
be solved in order to determine the unknown current distributions. The 
direct integration technique as outlined earlier is used to effect a nu- 
merical (discrete) solution for the current distributions. 

iicwreinstssitep in thes direct intepration technique is the choice of 
the expansion functions and the testing functions. As pointed out in an 
earlier section, the method of subsections in conjunction with point 
matching is a very convenient approach. The testing functions are Dirac 
delta functions since point matching is used. 

The expansion functions for the unknown currents are chosen to be 
piecewise constant for its ease of use and piecewise sinusoids for their 


improvement in convergence and required computation time. Symbolically, 
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the piecewise constant function is represented by 


N 
cap) = ie = gil) > WG = 


)] (58a) 


| sea S | ' ' | 
I (s ) - J 2IUG - Say - U(s' - a a (585) 


and where U(a) is the familiar unit step function 
U(a) = (59) 


Consistent with the above description, the current can be written as 


N-1 
t) — 7% ' \ zZ ' ' ' ' Z ' ' 
Itz ) = [U(z, =—7 + Lt, -Z") - U(z mre + [Uz ~Zy) (60a) 
and 
I (s') = 1“ws! om) ay 3p '_s!') U(s'-s' ee Ss.) 60b 
me )-1 (s,, - §s ANUXe 7): - U(s nae | ql BG all ) 


q=2 


Similarly, one can write for the piecewise sinusoids the following ex- 


pressions: 


Z Z, 
Nei sin k@* .-z') + I sin k(z'-z') 
T(z") = ) = nie te {U(z'-z") - WC sre 
n=1 sin K(z 7472.) c. 
n n 
(61a) 
S S 
Optiecan kis’ .—s") +° 1 sin k(s'-s') 
I (s') = ) q Gata qa q {U(s'~s') _ U(s'-s' 4) ] 
: q=l sin k(s' _-s') q q 
qtl q 


(GLb IE 


With the assumption that the segments were all equally spaced, the term 


2 is a constant. This allows a term delta (A) to be defined as 


1 
2 ———E——— 62 
: sin k (constant) ° see) 
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With this definition and equations (61) the current can be written as 





Z wh 
aL if 
a i to! 1_»! 1 : ia Soins 
I (2) A Sin K(z,-2 )U(Z, z') + sin k(z z Utz zy) + (63) 
N 
FL Z : ' ! fost = t,t : ! at tet t,t! 
A 2 Tatsin k(z a2) Ue Zi-1) Bes z'))+ sin K(Z0 4 Z ) (UCz zi) -Utz Zep) aI 


and similarly for the s current. With these expansion functions, the cur- 


rent distribution appeared graphically as shown in Figures 3 - 9, 
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Proeure 3. Real and imaginary parts of the current distributions for the 
radius = .001588\, theta = 0°, 36 segments. 
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Fagure 5. Real and imaginary parts of the current distributions for the 
Peis .0015ee,. theta = 0°, 30°, 45°, 60°, 85°, 90°, 36 segments. 


44 











AT O DEGREES 


X -KING 


—_ 
KK 
— 
x 
- 
60 
a 
vo 
ce 
° 
e 
uN 
N 
SITIAWIT [TTIW 
wa 
a 
G 
al 
a 
9 0 a 
>. =_ es, 
| oe ee Z 
= <£ O 
t [@) 
x ; 
fc 
x na 
= ~“ 
L | 
oo S 
ax 
=. 
ae) 
vo 
~~ 
x 
~ 
80 
i= 
ag 
wi 
O 
ron 
wm 
f 
iJ 
8 7, 
—=-— 
me 
[* 
es oO °o °o So ° o °o Lo) °o °o °o °o ° 
s . . e es e . se e* 4 e s 
Pee oe ay i rN “I = ° - “ " + wr 
1 t n ' t 1 
esatladmey{ it tw 


Figure 6. Real and imaginary parts of the current distribution for the 
radius = .003175A, theta = 0°, 30°, 45°, 60°, 85°, 90°, 32 segments. 
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Figure 7. Real and imaginary parts of the current distribution for the 
tadius = .004763\, theta = 0°, 20°, 45°, 60°, 85°,90°, 25 segments. 
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THETA = 45 DEGREES 
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Figure 8, Real and imaginary parts of the current distributions for the 
Bacbus = 001558), .0031/51, .004763i1, theta = 45°. 
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THETA = 90 DEGREES 
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Figure 9. Real and imaginary parts of the current distributions for the 
radius = .001588\, .003175\, .004763A, theta = 90°, 
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Vit RESUETS 


The system of integral equations and the proposed solution techniques 
are quite general, andthey allow variations in many of the parameters in 
the basic configuration. However, certain parameter variations appear to 
be more interesting than others. This study concentrates its efforts in 
three main areas: 

1. The ability of the numerical solution to give the result obtained 
by King in Ref. 26 when the angle between the two wires is zero. 

2. Current as a function of the angle between the two wires. 

3. Current as a function of wire diameter. 

The antenna is assumed to be driven by a delta gap generator. The 
current distributions are presented in terms of a dipole vice a monopole. 
The real and imaginary parts of the current are plotted for various bends 
of the antenna. The total length is \/4 for all the figures with the an- 
tenna bent at A/8 to give two equal lengths. Although these results are 
the only ones presented, the computer program allows the antenna to be 
bent at any point and the antenna to have any length. 

Figures 3 and 4 show the current for two different radii at an angle 
of zero degrees. for a comparison against King's frill feed in Ref. 26. 
These figures also show the result from Ref. 27 and 28 against King. 

References 27 by D.C. Kuo and B.J. Strait at Syracuse University uses 
a modified Pockington formulation with thin-wire approximations. 4 1s 
used instead of A and derivatives are approximated by finite differences 
over the same intervals used for integration. Galerkin's procedure is 


used with triangle expansion functions which are non-zero only over four 
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consecutive segments. Successive triangles overlap every two segments 
except on the ends of open wires. Junctions are treated as open wires 
each having a two segment overlap at the junction. An extended operator 
is used because the first order derivatives of the triangle functions 
feenot exist at certain points. 

Reference 28 by MB Associates and called the Antenna Modelling Pro- 
gram uses a Pockington formulation with thin-wire approximations. Point 
matching is used with the current being expanded as a sine, cosine and 
constant with an interpolation procedure so only N linear equations need 
to be solved. The extrapolated current for a given segment is forced to 
match the center current values in two adjacent segments. This interpol- 
ation procedure is also applied to wire junctions. 

Figures 5,6 and 7 give the currents for antennas of three different 
radii having 0, 30, 45, 60 and 90 degree bends. King's input values for 
a 0 degree bend are also plotted. Reece 8 and 9 show the currents for 
different radii antennas plotted on the same graph for bends of 45 and 90 
degrees, respectively. 

several general comments concerning the data presented may prove help- 
ful to the reader. All the results presented here are for the case of a 
one volt drive and assumes an open-ended perfectly conducting hollow 
cylindrical tube as the antenna. Extensive calculations have been carried 
out for bends of 0 to 90 degrees for each antenna radius given, but the 
plotting of all these angles on a single graph would obscure the overall 
view. It is believed that the data given is sufficient to give the reader 
an indication of the general variation of the current with respect to the 


parameter of interest, 
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VII. CONCLUSIONS 


In this thesis a new method is presented for analysis of an arbitrarily 
bent monopole. The arbitrarily bent monopole is a good choice because it 
contains cross—- coupling terms in the Hallén formulation which are not 
present in the O or 90 degree case. The currents had to be obtained 
using the appropriate boundary value problem based on Maxwell's equations 
which are formulated as a operator equation. The general procedure, 
"method of moments", is used to represent the operator by a matrix. Geo- 
metrical symmetries are utilized to simplify the matrix inversion. 

Numerical results are compared with data obtained by different theo- 
retical techniques and with experimental results where possible. Agree- 
ment is generally very good and this is interpreted as an indication of 
the accuracy of the solution. 

Certain parameters should now be concentrated onif the conclusions 
are to be used properly. These efforts should be in these areas: 

1. Choice of bases 
2. Antenna drive 
Se «6Pascontinuity at the bend. 

As can be seen by equation (36) the sinusoid current expansion 
equations are continuous but their derivatives are not. This implies 
mathematically that there are pockets of charge at the end of each 
subinterval which physically are not there. If one is interested in 
the near field, these mathematical charges cause discontinuities in the 
near field. One solution has been to use a least square fit to the cur- 


rent values obtained. 
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The delta gap drive was chosen because it was an ideal model. Even 
though the delta gap gives a physically meaningless infinite input sus- 
ceptance associated with the knife edges, the answer is close to that 
achieved by experiment. The reason is that any method that approximates 
the current by a series of continuous functions will not reproduce the 
extremely sharp rise in amplitude near z = z' unless a very large number 
of terms is used. A better drive would be the frill feed as given in 
reference [29] which corresponds closely to the actual drive used on 
physical antennas. 

Although not plotted on the figures given, a discontinuity in cur- 
rent was obtained within one subinterval either side of the bend. This 
was to be somewhat expected because the formulation considers two straight 
wires joined at a point which gives a sharp bend where the cross section 
is abruptly changed. When one uses the approximate kernel, I(z) must be 
a slowly varying function in the vicinity of z = z' for the approximate 
kernel to be within an order of the square of the radius of the exact 
kernel. This is not true at the bend. It is this slow variation 
or fast variation that causes the various methods of iteration (Hallén, 


Gray, King-Middleton) to require more or fewer terms in representing the 


mua rent. 
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APPENDIX A: REDUCED KERNEL 


Equation (39) in the thesis must be investigated and R and P deter- 
mined in a more precise manner. A detailed development of the vector 
potential must now be presented. 

The vector Sonentlan, ate to a volume current density can be written 
as 
~jk[r-r'| 


dV(r') 


n@y aan {| J(x'Je (1) 


Ant fo aa 

yt 16 ie 
Looking at Figure 10, it is seen that the vector r locates the point at 
which A(r) is to be determined and r' locates the source point. With a 


symmetrically driven antenna whose wall thickness approaches zero, the 


vector potential has only a z component which can be written as 


— ul o ° Es Gay nye Jkie-r' | 
hee) Eas | | | SA oe dy do'daz' (2) 
Z +77 lr 7 z| 
z' ¢'=0 r'=0 
or integrating over r' gives 
Zt — — 
~jk|r-r'| 
aa tl) ValieiGar 2) lie ee 
ANG) Sie | | = (3) 
z' o'=0 (es mae 5 
Now due to symmetry 
Ca) 2taJ (a,2") (4a) 
OF ! = alk 1 
aJ (a,z Se Gy. (4b) 
Substituting (4b) into equation (3) gives 
u il “Tg [Foz | 
= —- ' pan eee | 1 1 5 
eal | Me Nig all pe 2 
Za o'=0 


ae 








Figure 10. Antenna Geometry assuming a surface current in 
tales GLreCtiOn. 
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where r' is constrained to remain on the surface of the cylinder. 
Appealing to basic vector algebra and trigonometry and being aided 


by Figure 10 one can write r and r' in the following form: 
eel ' cee [ ah t ug [ 
rs x, a vil + ZU, mecos ou, Casco UE ZU, (6a) 
Pee eee z= tr eos ou, - r' sin oe ate ae COD, 


Now |x -r'| is written as 
| 1 
fe Feo], Sa Care Rb (y = yt)? + (= 2')*)? 
il 
ee + Ca - 2rr'cos(¢ - ¢') + (z - Beye (7) 


I 


Lf |r - r'| is evaluated at r' equal to the radius of the z-directed 


wire and @ equal to zero then 


Ik 
ie sai ae = ie + ae - 2ra cos ¢' + (z - oe (8) 
g =0 
Trigonometry gives 
2 
1 = cos 28 = 2 sin 8B (9) 
Applying equation (9) to (8) with r = a gives 
i 
— —_ q 
| x - a} - eg sin” (%) + (z - Bene : (10) 


Going back to equation (5) and inserting equation (10) gives 
2 ,¢' 2 
: 1 favak[4a’sin’ 2? + (2 - 2!) fee 
Pal Lae { aoe DES ea Aes, a cae ea J 1 
mn i | eer oats eo + 222 acl a 
ot 4"=0 a sin G Z 
One sees immediately that at z equal z' and-¢'equal.0O or 27 the inner 
integrand is singular; however, the integrand is logarithmically singular 
and is integrable. Even though it is integrable, evaluation of the in- 


tegrand by numerical techniques is extremely difficult. This is one 


major reason why a reduced, average or approximate |r - ie” is found. 


a2 





For |z - z'| >> a and for thin wires where a << \ the following 


approximation can be obtained. 


2/2 yg 


z= 7] = [2 - a") + @a sin (*)) =< (12) 


Even this for thin wires is not a good approximation where z = z'. Ob- 


| 

serving that 2a sin e) Vemtheeaistance from a point (a,0,z') to a 
point (a,¢',z), it is seen that both points are on the same circle of 
radius a in a constant z' plane. Looking at Figure 10 it is seen that 

d' , 
Zee sin ) was a distance that assumed values from 0 to 2a as 9¢' 
varied from 0 to 2m with a mean value of a. Placing the mean value into 
equation (11) and integrating gives 
2 ele 


e eee + (z - z') 


A (a,z) = = | I(z')= (13) 
: su oi faa a (zZ = z") 


1 
2,1/2 oe 


It is restated that equation (13) is valid only for a << i and 
a << h. Now going back and using Figure 11, equation (5), and assuming 
the current to be a filament along the z axis, the vector potential 


A (r) at a point (a,0,z) would be 


° _-akla® ee CZ = Beg 

— a Ne : 

A_(a,0,2) Te (ez) 5 POLE dz (14) 
ae fay =-ec2)—) 2) 

which is exactly the same equation as equation (13). This is a most 


important point. Using the reduced kernel is the same as treating sur- 
face current on the cylinder as if it were an equivalent line current 


along the cylinder axis. 
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Figure 11. Antenna Geometry assuming a filamentary current 
along the Z axis. 
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APPENDIX B: SOLUTION OF COUPLED INTEGRAL EQUATIONS 


Equation (10) and (11) in the thesis give the two coupled partial 
differential equation that must be solved. Equation (12) gave the form 


that equation (10) and (11) must be in to apply the solution given in 





equation (17). Equation (10) gives the form for the z-directed wire as 
52 _ 12 
74, @] + A) = aha) - Sta) | (10) 


Notice that A (x) in cylindrical coordinates depends on (R,$,z) where 
R is the radius vector in these coordinates. If one fixes R and ¢, then 
A (x) will be a function of z alone along any line parallel to the z-axis. 


Equation (10) then becomes 


2 KP : 
ae Gye A (2) = (2) = sala A, (2) : (1) 
Pawation (12) in the thesis is 
Ove 
ye. Ho E(x) = U(x) = 22 (12) 


ax 


Examining equation (12) it is seen that equation (1) may be substituted 
directly into equation (12) and equation (17) which is the solution to 
equation (12) may be applied to equation (1). Thus it is seen that for 
equation (10) to use the solution given by equation (17) it must be in 
equation (1) form. This means equation (10) Palhition is good only along 


a line parallel to z-axis. 
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APPENDIX C: REDUCTION OF A DOUBLE INTEGRAL 


In order to show the solution to equation (47) in general, the sol- 
ution will be merely given and then differentiated to show that it does 
indeed give equation (47) even with arbitrary signs and the addition of 


mmeconsctant. the solution is 


: = -jkR 
{ JE + Ze ee k(B-e) de eee (= cos (ee) - jRsin k(6-e) } 

R R R[R-(Ce-f) ] 
where (1) 

Z Zz Pevaley 7, 
Peemia + 6 (c) + fe-£(<)) } i 
Now let the solution in general be 
27 JKR 
= { (e+f£+c) cos k(8-e). - jRsin k(B-e) } (2) 
R(R-[e+f(n)tec] 7) 
and 
R = les + ey + [exE(n)te]-}o!“, (3) 
where g and f are functions of the parameter n. 
Now taking the derivative of the function 
q a JKR 

famous | (et Ftc) cos k(B-e€) - jR sin k(6-e)]} = 

R -R(¢etftc) 

F -~4kR 
[(exftc)cos k(B-e) - jR sin k(B-e)] AS a See sh 
R -R(etftc) 
2 SE d 
ap a je [(etftc)cos k(B-e) - jR sin k(8-e)]. (4) 
R -R(e+tftc) 
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Working on the first derivative in equation (4) gives 


~jkR ~ a 
l=] = [ir -rterete) Hanele 
R -R(etftc) R 
—{kR 3° ( rc) (sari tas ee: ) 1 
2 ES [ee SSSR | ONG) o See a genes, 
R R [R°-R(etftc) ] 


This reduces to give for the first part of equation (4) 


TR spec) GE - 
—=- 4°, |(etftc) cos k(f-e) - jR sin k(B-e)] . (>) 


Ro eee 


Similarly for the second derivative in equation (4) gives the following 


result 


[(etftc)k sin k(B-e) + cos k(B-e) + jkR cos K(B-e) - j_sin tire) teste) 


-4kR 
a ae (6) 


PEeR(e+ftc)— 


Adding equations (5) and (6) and factoring out the common elements re- 


sults in 


~jkR 
5 = [(Cetftc)k sin k(B-e) + cos k(B-e) + jkR cos k(B-e) - (7) 
R —-R(etftc) 
2 ss il 
j(etftc)sin k(B-e) x (etftc) (-jk-R)cos k (Bre), (etftc) (-jk-R) (-jR sin A 
R R R 


Now start cancelling common terms 


Cetftc)k sin k(f-e) - (etftc)k sin k(g-c) = 0 (8a) 
Weete eam Kk (6-2) ii(etiic)sin k(B-c) _ 
| R cies 0 (8b) 
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which gives as a result after some algebra 


Ags cos k(f—-e) 


ee E. Mee=(e2esey Til + qkRI}) (9) 
R -R(etfitc) R 


Cancelling common factor, one has 


eek vee af 
—__. 2 a + jkR)} =e apd ee (B-c) its + 1 (10) 


R R R R 
which is equation (47) in the thesis. 
Thus equation (2) and (3) are required functions in general and 
equation (1) is the necessary function when there is no arbitrary con- 


Stant or arbitrary signs in the integrand. 
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